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We present a theoretical survey of ferromagnetic transition temperatures in cubic (III,Mn)V semi- 
conductors based on a model with S — 5/2 local moments exchange-coupled to itinerant holes in 
the host semiconductor valence band. Starting from the simplest mean-field theory of this model, 
we estimate the Tc enhancement due to exchange and correlation in the itinerant-hole system, and 
the Tc suppression due to collective fluctuations of the ordered moments. We show that high critical 
temperatures in these ferromagnetic semiconductors require both the large magnetic susceptibility 
contribution from the valence band's heavy holes and the large spin stiffness that results from strong 
spin-orbit coupling in the valence band. Our calculations demonstrate that this model for the ferro- 
magnetism of these systems is fully consistent with the room-temperature ferromagnetism reported 
for Mn doped nitrides. 



The 1992 discovery^ of hole-mediated_ferromagnetic 
order in (In,Mn)As has motivated researcho on Mn doped 
GaAs and other III-V_host materials. Ferromagnetic 
transition temperaturesa in excess of 100 K and long spin- 
coherence time^ in GaAs have fueled hopes that a new 
magnetic medium is emerging that could open new path- 
ways for information pijocessing and storage technologies. 
Recently, confirmatiatfl of the room-temperature ferro- 
magnetism predictedQffl in (Ga,Mn)N has added to inter- 
est in this class of materials. 

When substituted on the cation site of a III-V semi- 
conductor, it is generally accepted that Mn acts as an 
acceptor leaving a Mn^+ half-filled d-shell with angular 
momentum L — and spin S = 5/2. There is also exper- 
imental evidenceErtl that ferromagnetism occurs in these 
materials because of interactions between Mn local mo- 
ments that are mediated by holes in the semiconductor 
valence band. The local moments and the itinerant holes 
are coupled in this approach by an antiferromagnetic ex- 
change interaction. The participation of itinerant holes 
in the ferromagnetism of these diluted magnetic semicon- 
ductors adds to their richness, leading to strong magneto- 
transport effects that might have important applications 
and, of particular interest to us here, to transition tem- 
peratures that are sensitive to the details of the host band 
structure. 

The description of ordered states in (lll,Mn)V systems 
is greatly simplified by a virtual-crystal type approxima- 
tion in which the Mn ion distribution is reuLaced, by a 
continuum with the same impurity density.Hil3lij Dis- 
order effects associated with randomness in the Mn ion 
sites can then be treated perturbatively, when necessary. 
(For example randomness in the Mn distribution likely 
limits the conductivity at low-temperatures.) This strat- 
egy will fail in the limit of dilute Mn ions and also when 
the exchange interaction between band and Mn spins is 



too strongEi It does, however, seem to be reliable in the 
limit of principal interest, that of high Mn densities and 
high critical temperatures, where the holes are metallic 
and their interaction with Mn acceptors will be effectively 
screened. 

There is at present considerable activity directed to- 
ward the growth of many different (III,V) materials con- 
taining Mn. Inlhis paper we present and discuss theoret- 
ical predictionsE3 for their ferromagnetic critical temper- 
atures based on the continuum model. We go beyond the 
standard mean-field theory of this model by accounting 
for the role of Coulomb interactions holes in the valence 
band, which enhances the critical temperature, and for 
correlations in Mn ion orientations which reduce the ener- 
getic cost of unaligned spin coE^suxations and therefore 
lower the critical temperature.EZrt^l 

When disorder and hole correlations are neglected, the 
meaxLjfield transition temperature of this model is given 
byBS 
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where A^Mn — 4a::/aj*j. is the Mn density in Mnajllli-^V 
zinc-blende semiconductors with a lattice constant aic, 
Jpd is the localized spin - itinerant spin exchange cou- 
pling constant, and x/ is the itinerant hole magnetic sus- 
ceptibility. To understand the qualitative implications 
of this Tc-equation (|^), we first consider the case of a 
model itinerant electron system with a single spin-split 
band and an effective mass m*. The kinetic energy of 
the band holes gives a contribution to the susceptibility: 
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m*kF 
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where kp is the Fermi wavevector. For weak interactions 
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the exchange energy of the spin-polarized parabolic-band 
model adds a contribution: 

where e is the dielectric constant of the host semicon- 
ductor. At high hole densities p, the kinetic-energy term 
dominates and the mean-field critical temperature, , 
is proportional to the Fermi wavevector, i.e., to p^^^. 
Equations and (||) also show that the band contri- 
bution to the susceptibility increases linearly with m* 
while the exchange correction is proportional to (m*)^. 
(Using this simple band model as a guide, hole correla- 
tion effects, which will not be discussed in detail here, 
suppress the critical temperature by several percent for 
typical experimental hole densities, p ~ 0.1 nm~^.) 

To obtain quantitative predictions for the critical tem- 
perature, it is necessary to evaluate the kinetic and ex- 
change contribution to the itinerant hole suscppjihility 
using a realistic six-band Kohn-Luttinger modelErQiHiJ, in- 
stead of the parabolic-band model. The band Hamilto- 
nian contains the spin-orbit splitting parameter Ago and 
three other phenomenological parameters, 71, 72, and 73, 
whose values for the specific III-V host can be found, 
e.g., in Refs. 0,^. Only zinc-blende crystals are con- 
sidered here; in nitrides, which crystallize in both zinc- 
blende and wurtzite structures, the critical temperature 
is expected to bs-|0nly weakly dependent on the lat- 
tice configuration.Q The non-interacting-hole mean-field, 
, and the exchange-enhanced, T^^, critical tempera- 
tures for GaAs and GaN host semiconductors doped with 
5% of Mn are plotted in Fig. |l| as a function of the hole 
density. Transition temperatures for other III-V hosts 
and for hole densities p = 0.1 and 0.5 nm~^ are listed 
in Table ||. We,|have assumed Jpd = 55 meV nm"'^ for 
all III-V hosts.ElEil In the density range considered, only 
the two heavy-hole and two light-hole bands are occu- 
pied in arsenides and antimonides. However, the mixing 
between these four bands and the two spin-orbit split-off 
bands is strong and must be accounted for. In nitrides 
and phosphides, spin-orbit coupling is weaker and all six 
bands are occupied by holes. The numerical results for 
the transition temperature are consistent with the quali- 
tative analysis based on the parabolic band model: 
follows roughly the p^^^ dependence, the ~ 10 — 50% 
exchange enhancement of the critical temperature has a 
weaker density dependence. j— , 

For the thoroughly studied (Ga,Mn)As matcrialEll 
doped with 5% of Mn and with p = 0.35 nm^'^, we 
find T^^^ = 102 K and T^"" = 112 K. This result is 
in remarkably good agreement with the experimental 
transition tempetature of 110 K . Good agreement be- 
tween experimenttll and mean-field theory is also found in 
(In,Mn)As compounds, that have lower critical temper- 
atures compared to their (Ga,Mn)As counterparts, pri- 
marily because they have smaller band masses. An es- 
pecially important system is (Ga,Mn)N, which appears 



to exhibit room-temperature ferromagnetismU Our cal- 
culations for this III-V semiconductor give « 400- 
600 K and T^"" ^ 600-900 K for Mn content x = 5% over 
the range of hole densities studied. Although we cannot 
compare those numbers directly to the igcently measured 
Tc = 228-370 K in (Ga,Mn)N samplesl, since they have 
inhomogeneous doping profiles (resulting from the chem- 
ical vapor deposition growth technique), the mean- field 
theory predictions have the correct order of magnitude 
but appear to overestimate the transition temperature. 

We now establish the quantitative reliability of the 
mean-field theory in arsenides (and also in antimonides 
and lower-carrier-density phosphides) and explain why 
mean-field theory overestimates the transition tempera- 
ture in the heavy effective mass nitrides. Our arguments 
are based on the estimates of the suppression of the tran- 
sition temperature due to collective spin-wave fluctua- 
tions. We emphasize below that mean-field theory would 
fail badly in any of these materials if the appropriate 
band model were really a single parabolic band with a 
mass equal to the material's heavy-hole mass. Strong 
spin-orbit coupling in the valence band, which yields 
band states whose orbital and spin character changes 
completely|-as the band wave-vector varies over the Bril- 
louin zone£3 enhances the spin stiffness, increases the en- 
ergy of correlated reduced magnetization configurations, 
and reduces the imcortance of corrections to the mean- 
ficld-theory result .Ej 

Isotropic ferromagnets have spin-wave Goldstone col- 
lective modes whose energies vanish at long wavelengths, 

rik = Dk^ + 0{k^) , (4) 

where k is the wavevector of the mode. Spin-orbit cou- 
pling breaks rotational symmetry and leads to a finite 
gap. According to our numerical studies J13 this gap is 
small however, much smaller than Tc for example, and 
plays a role in magnetic fluctuations only at very low tem- 
peratures. Spin-wave excitations reduce the total spin 
by one, at an energy cost that is, at least at long wave- 
lengths, much smaller than the mean-fleld value, JpdS^"^ 
where s'-^^' is the mean-field band spin-density. The im- 
portance of these correlated spin excitations, neglected 
by mean-field theory, can be judged by evaluating an 
approximate Tc bound based on the following argument 
which uses a Debye-like model for the magnetic excitation 
spectrum. When spin-wave interactions are neglected, 
the magnetization vanishes at the temperature where the 
number of excited spin waves equals the total spin of the 
ground state: 

NMnS=^ J^\kkM^k), (5) 

where n{Qk) is the Bose occupation number and the De- 
bye cutoff, k]j — {Gn'^NMny^^- It follows that the ferro- 
magnetic transition temperature cannot exceed 

AbT™" = ^^klD{T^°'') . (6) 
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In applying this formula to estimate Tc we have approx- 
imated the temperature dependence of the spin stiffness 
by 

D{T)=Do{S){T)/S , (7) 

as we did in ReL-. |23[ where Dq is the zero- 
temperature-atiffnessjlj and {S){T) is the mean-field Mn 
polarizations at a temperature T. If the difference be- 
tween T™'' and is large, the typical local valence- 
band carrier polarization will remain finite above the crit- 
ical temperature and ferromagnetism will disappear only 
because of the loss of long-range spatial order, the usual 
circumstance for transition metal ferromagnetism for ex- 
ample. 

In discussing corrections to mean-field-theory Tc esti- 
mates, we compare spin-stiffness results obtained with 
the simple two-band and realistic six-band models. De- 
tails on the formalism used to calculate Dq can be found 
in Refs. We find that the zero-temperature spin 

stiffness is always much larger in the six-band model. For 
(Ga,Mn)As, for example, the two-band model underesti- 
mates Dq by a factor of ~10-30 over the range of hole den- 
sities considered. Furthermore, the trend is different: in 
the two-band model the stiffness decreases with increas- 
ing density, while for the six-band description the initial 
is followed by a saturation. Even in the limit of low car- 
rier concentrations, it is not only the (heavy-hole) mass 
of the lowest band which is important for the spin stiff- 
ness. In the realistic band model, heavy-holes have their 
spin and orbital angular momenta aligned approximately 
along the direction of the Bloch wavevector. Exchange 
interactions with Mn spins mix the heavy holes with more 
dispersive light holes. Our calculations show that heavy- 
light mixing is responsible for the relatively large spin 
stiffnesses which are responsible for the general success 
of mean-field theory. Crudely, the large mass heavy-hole 
band dominates the spin susceptibility and enables local 
magnetic order at high temperatures, while the dispersive 
light-hole band dominates the spin stiffness and enables 
long-range magnetic order. The multi-band character of 
the semiconductor valence band plays an essential role in 
the ferromagnetism of these materials. 

In the insets of Figs. |^(a) and (b), solid lines show 
the spin stiffness for GaAs and GaN, doped with 5% of 
Mn and with p — 0.4 nm~^, as a function of temper- 
ature. The critical temperatures, T™'', were obtained 
from these curves and from Eqs. d) and @). T™", cal- 
culated using a non-interacting hole approximation, and 
r°'^, calculated including exchange interactions but based 
on a mean-field treatment of the Mn spin-orientation 
thermodynamics, represent approximate lower and up- 
per bound for the ferromagnetic transition temperature 
in (III,Mn)V semiconductors. With the exception of ni- 
trides and higher- carrier- density phosphides, exchange 
enhancement and correlated fluctuations change Tc by 
less than 20%, and these two influences will tend to can- 
cel. This property explains the remarkable success of the 



mean-field theory in these itinerant ferromagnets. In the 
heavy effective mass nitrides, mean-field theory looses its 
quantitative accuracy, but still predicts the correct order 
of magnitude for the ferromagnetic transition tempera- 
ture. 

The dot-dashed lines in Fig. |^, and Table | summa- 
rize critical temperature estimates T^**' obtained from 
Eqs. (H) and (Q) and by shifting the _D(T)-curves in the 
insets of Fig. ]l[ along x-axis by T^"" - T^^. This ap- 
proach accounts, in an approximate way, for both ex- 
change interactions and correlated spin fluctuations. For 
arsenides, antimonides, and lower-carrier-density phos- 
phides, Tc^^ is very close to T^^ justifying the mean- 
field description of ferromagnetism in these semiconduc- 
tor compounds. The estimated critical temperatures for 
(Ga,Mn)N, = 200 - 400 K, refiect a more significant 
suppression of ferromagnetism due to collective excita- 
tions. Our theoretical T"^^ are in good agreement with 
available experimental transition temperatures in stud- 
ied (III,Mn)V diluted magnetic semiconductors, includ- 
ing (Ga,Mn)N. 
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TABLE I. Mean-field {T^^), exchange enhanced (T^), 
collective (T^°"), and estimated (T^^*) ferromagnetic transi- 
tion temperatures in III-V host semiconductors doped with 
5% of Mn and with itinerant hole densities p = 0.1 and 



p — 0.5 nm 
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FIG. 1. Main graphs: Mean-field (T^^), exchange en- 
hanced (TD, collective (T™"), and estimated (T^*) ferro- 
magnetic transition temperatures in GaAs (a) and GaN (b) 
host semiconductors doped with 5% of Mn are plotted as a 
function of itinerant-hole density. Insets: temperature de- 
pendence of the spin stiffness (solid lines) calculated from the 
zero-temperature stiffness and from the mean-field magneti- 
zations. The hole density p — 0.4 nm~^ and x = 5% for 
both GaAs (a) and GaN (b) hosts. Shifted (dotted) curves 
account for local magnetization enhancement due to exchange 
interaction in the itinerant hole system. 
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